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Abstract 

We propose a new real-space renormalization group transformation for dy- 
namical triangulations. It is shown to preserve geometrical exponents such as 
the string susceptibility and Hausdorff dimension. We furthermore show evi- 
dence for a fixed point structure both in pure gravity and gravity coupled to 
a critical Ising system. In the latter case we are able to extract estimates for 
the gravitationally dressed exponents which agree to within 2 — 3% of the KPZ 
formula. 



1 Introduction 



The success of dynamical triangulations (DT) in providing a discrete regularisation 
for two dimensional quantum gravity is well known [0. The advantages of this 
approach are twofold; firstly it allows powerful analytical techniques based on matrix 
models to be employed and secondly it lends itself to non-perturbative studies via 
numerical simulation. The latter approach can be applied both to study models 
which are currently not amenable to exact treatments (for example with matter 
central charge c > 1) and to compute quantities of physical interest, such as the 
nature of the quantum geometry, which are difficult or impossible to arrive at with 
analytical methods. 

Much of the numerical work has centered around the application of finite 
size scaling to extract critical exponents (see e.g. 0). In statistical mechanics the 
validity of such techniques depends on the existence of a renormalization group 
(RG) which governs the flow of couplings under changes in scale. It is tempting to 
believe that a similar RG structure must underlie the DT models and is responsible 
for the observed scaling. In the continuum formulations of quantum gravity the 
very issue of a renormalization group is a difficult one to formulate since the theory 
(in the absence of a cosmological constant) possesses no length scale. In contrast 
the DT formulation contains an invariant cut-off corresponding to the elementary 
triangle edge length. The latter may be traded in for the number of triangles N if 
the physical volume is held fixed. 

A successful RG transformation would be important both conceptually and 
as a powerful new tool with which to compute critical points and critical exponents 
for systems coupled to quantum gravity. Recently various real-space RG transfor- 
mations for dynamical triangulations have been proposed. The first of these ap- 
proaches, due to Renken achieves this by randomly selecting a subset of nodes 
on a given triangulation. These nodes are then connected by links to form coarser 
triangulations, in such way as to preserve locally the relative geodesic distances of 
the selected nodes. We can call this a local geodesic blocking (LGB). 

Another approach advocated by Krzywicki et al. Q constructs a coarser tri- 
angulation by eliminating a class of extremal baby universes associated with the 
original random lattice (extremal implies that no further baby universes grow on 
them). Since the baby universe distribution is intimately connected to the fractal 
structure of the ensemble of triangulated manifolds, it is hoped that this transfor- 
mation will (approximately) preserve the fractal structure. This may be termed a 
fractal blocking (FB). 

In this paper we propose a new real-space RG transformation for dynamical 
triangulations based on a local node decimation scheme (ND). The essence of this 
method is that nodes are removed randomly from a triangulation. After each one is 
removed its neighbors are reconnected in such a way as to preserve the integrated 
"local" curvature and maintain the triangulated structure of the manifold. As we 
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will demonstrate this method preserves the relevant long distance physics of the 
surfaces, and in addition allows us to compute the critical behavior of matter fields 
living on the ensemble of random manifolds. 

The paper is organized as follows: In section 2 we introduce the methodology 
of Monte Carlo renormalization group. In section 3 wc discTiss the problems involved 
in blocking dynamical triangulations and describe the method we propose. In section 
4 we apply the method to pure gravity, demonstrate that it preserves the fractal 
structure and show evidence for a fixed point behavior. We also look at the effects of 
including irrelevant operators in the action. In section 5 wc apply the method to a 
critical Ising model coupled to gravity and show that, combined with an appropriate 
blocking of the Ising spins, it yields the correct critical behavior. In section 6 we 
discuss some alternative blocking schemes. Finally section 7 summarizes our results. 

2 Monte Carlo Renormalization Group 

In conventional lattice field theory a RG or block transformation acts so as to 
reduce the number of degrees of freedom by replacing the original (bare) theory by 
one defined on a coarser lattice. This effectively integrates out the short distance 
fluctuations, allowing us to focus directly on the long distance physics which governs 
the critical behavior. Associated with this blocking there will be a flow in the 
coupling constants of the system to compensate for the change in length scale. If 
the system is close to a continuous phase transition this flow will pass close to a 
critical fixed point. Studying the flow in the vicinity of such a critical point allows 
us to determine the critical exponents of the model. 

Consider a system of spins {a} on the sites of some regular lattice with 
interactions described by a Hamiltonian H(a). All thermodynamic quantities can 
be found from a detailed knowledge of the partition function 

Z{H) = Y^eHi^). (1) 
{<^} 

We define a RG transformation that maps the system onto a blocked system, with 
fewer spins ji, by 

e^'(^) =5^ P((7,/x)e^(^). (2) 

{^} 

P(c7, ji) is a projection operator that couples together the spins of the original and 
the blocked system. Note that as the background lattice is regular, the coarser 
lattice trivially has the same structure. 

If the RG transformation is to preserve the physics of the model the projec- 
tion operator should satisfies the following conditions; 

P(c7,At)>0 V (T,/x and Y^P{a,n) = l. (3) 
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The former condition ensures that the renormahzed Hamiltonian is real, the latter 
that the partition function is preserved. In addition it is usually required that 
P{(T, /i) preserves the symmetries of the Hamiltonian. 

Expanding the Hamiltonian on some operator basis H(a) = J2a -^aOa{cr) 
the definition of a RG transformation Eq. || implies a set of recursion relations for 
the couplings 

kC^+i) =i2[KW]. (4) 

The critical behavior of the model is calculated by looking at fixed points of these 
recursion relations, K* = R[K*]. Close to such a fixed point the flow of the dif- 
ference in the couplings from their fixed point values may by analyzed in a linear 
approximation; 

'^^^ - (5) 

K=K* ^ 

The critical exponents are then determined from the eigenvalue equation 

We distinguish between three types of eigenvalues: Aj > 1 is called a relevant eigen- 
value, Aj = 1 marginal and Aj < 1 irrelevant. The meaning is simple. Consider 
expanding the operators in the effective action on a basis given by the eigenvectors 
of Eq. (6). Then, under iteration of the RG step, components along the direction 
of an irrelevant eigenvector (Aj < 1) will be driven to zero. Conversely, components 
along the relevant directions will increase. These directions then describe how the 
renormahzed Hamiltonian flows away from the fixed point. The case of a marginal 
eigenvalue is special and has to be studied closely in any individual case. 

A relevant eigenvalue defines a scaling exponent yi = log A^ / log bi which 
describes the singular behavior of the free energy [bi = A^(^) /Ar(fe+i) ig the volume 
rescaling). It can be related to a combination of critical exponents of the model. 

Although the condition imposed on the projection operator Eq. (^) allows 
plenty of freedom in defining a RG transformation, in practice it is very hard to 
find one that allows for an explicit calculation of the linearized matrix Tq,^. One 
solution is to combine an exact RG transformation with Monte Carlo simulations 
- the Monte Carlo renormalization group. Monte Carlo simulations are used to 
generate a sequence of configurations characteristic of the original bare Hamiltonian, 
which then are blocked using the RG transformation. If the transformation is 'apt' 
this will result in a configuration describing the same long distance physics, but on 
a coarser lattice. There are two main sources of error in such a calculation; the 
statistical error arising from the finite number of configurations generated by the 
Monte Carlo algorithm and the necessary truncation of the operator basis for T^/s- 
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The former can be estimated and systematically reduced, the latter can be assessed 
by varying the operator basis. 

The renormalized Hamiltonian is approximated using a set of operators O-^ 
on the finite lattice and the linearized matrix Tap calculated by solving a set of chain 
rule equations: 

d < o\'^'^ > ^ d < oy+^^ > 



where 



d < oi') > 



(fc) ^ 

J-^ = < 0('=) 0(f) > - < 0(^) X 0(f) > . (9) 



The critical exponents are then calculated from the eigenvalues of T^/s- One ap- 
pealing feature of this method is that it is relatively easy to include a large set of 
operators, allowing for a systematic improvement of the approximation. 



3 Random surfaces 

The model we are interested in is the dynamical triangulation formulation of 2d 
quantum gravity. It is defined as a sum over all possible ways of gluing equilateral 
triangles together in such way that they form a closed (piecewise linear) manifold 
T; 

Zifi) = E e-'^^ZM(r). (10) 
TeT 

IX is the cosmological constant, the volume of the triangulation and T a suitable 
class of triangulations. In this paper we chose this class to be combinatorial two- 
manifolds, i.e. no two nodes are connected more than twice and no node is connected 
to itself. This implies that there are no nodes of order two or less. Finally Zm (T) 
denotes the partition function corresponding to a sum over all possible matter fields 
living on the triangulation T. 

In numerical simulation it is more convenient to study this model at fixed vol- 
ume, i.e. to look at the micro-canonical partition function Z{N) related to the grand 
canonical partition function Eq. ( p!o|) through a discrete Laplace transformation 

Z(^) = Ee-^^Z(iV). (11) 
JV 

To construct a blocking transformation for this model is quite different than 
for usual field theory models on regular (flat space) lattices. There the RG trans- 
formations trivially preserve the lattice structure and the main problem is to devise 
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reconnect 




Figure 1: The five possible ways of re-triangulate after removing a node of order 
five. Choosing one of these possibilities at random yields a blocked triangulation. 
If a node of order six is removed there are 42 ways of reconnecting. 

a blocking scheme for the matter fields. In contrast it is not possible to preserve ex- 
actly the features of any given fine triangulation under blocking. Indeed the lattice 
itself now carries dynamical degrees of freedom. Thus any transformation will have 
to replace a given triangulation T{N) by some triangulation T'{N'), corresponding 
to a blocking factor bi = N/N' , 

T'{N') = R[T{N)]. (12) 

Clearly the choice of the transformation R is of crucial importance. Presumably it 
should have the property of preserving certain aspects of the long distance geom- 
etry. But what aspects? We propose a very simple method based on a local node 
decimation of the triangulation. 

First pick a node at random in the initial triangulation T. If we attempt to 
remove this node then unless its coordination is three we will be left with a polygonal 
"hole" in the triangulation. However by suitable link addition it is possible to re- 
triangulate the interior of this hole. This can be done in many ways; we choose one 
of those at random. This is illustrated in Fig. ||. 

By iterating this procedure an arbitrary number of times we can produce 
a blocked triangulation of any volume. In practice, this procedure is effected by 
randomly flipping the links around the selected node until its coordination number 
becomes equal to three. At this point any curvature associated with the node has 
been smeared out over its neighbors. We then remove the three-fold node. It is 
trivial to see that this method preserves the integrated curvature within a sufficiently 
large loop around the node in question. However it is not at all obvious that this is 
enough to capture the most important long distance physics governing the fractal 
structure, critical exponents etc. We examine these issues now. 
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4 Pure gravity 



We have applied this blocking method to the case of pure gravity (triangulations 
taken with weight unity) and performed a series of Monte Carlo simulations on 
lattices ranging from 250 to 4000 nodes. A standard link flip algorithm was used 
to explore the space of configurations and typically 10^ sweeps performed for each 
lattice volume. Every twenty sweeps a configuration was picked out and blocked 
several times down to 8 nodes. 



4.1 The fractal structure 

To see whether the long distance physics of the surfaces is being preserved we have 
measured two exponents which describe complementary features of the fractal struc- 
ture; the string susceptibility exponent 7^ and the HausdorjJ or fractal dimension dfj- 
The string susceptibility exponent is related to the large volume behavior of 
the fixed area partition function Z{N), 

Z{N) ~ N^^-^ef"^^. (13) 

In the case of pure gravity 7^ = —1/2. An efficient method for measuring 7^ in 
numerical simulations is to look at the distribution of so-called minbu's (minimal 
neck baby universes) ^ . A minbu is a part of a triangulation which is connected 
to the rest through a neck or a loop consisting of three link^. The distribution of 
minbu's can be calculated from the number of ways a baby universe of volume B 
can be glued onto a surface of volume N — B; 

It is easy to measure this distribution numerically and the exponent 7^ is obtained 
by fitting to Eq. (p!4|). The accuracy of the fit can be improved considerably by 
including a 1/-B correction term. 

We have measured 7^ in this way for the ensemble of triangulations T^^^ 
obtained by blocking k times. The result is shown in Fig. ^ where we plot 7^^ 
versus volume. This is done for four different initial volumes A^^''). It should be 
noted that the number of blocking steps increases to the left in the plot. It is clear 
from Fig. || that independent of the initial volume and number of blocking steps we 
get the correct value of 7^ down to volumes of size 200 nodes. For smaller volumes 
it is hard to measure 7s accurately due to finite size effects. This demonstrates that 
even after removing 90% of the initial triangulation by blocking the model remains 
in the universality class of pure gravity. 

Additional evidence is obtained by measuring the Hausdorff dimension of the 
surfaces. It has recently been shown that the Hausdorff dimension is related to point- 
point distributions on the triangulations and that it can be measured numerically 

^As we do not allow degenerate triangulations this is a minimal loop on the surface. 



6 



(k) -0.35 



-0.4 



-0.45 



-0.5 



-0.55 



[] 



o 



100 



[] 
)( 



h =^2 



o ^v^'^^ 


= 250 




= 500 




= 1000 


X 


= 2000 



Pure gravity 



1000 



N 



(k) 



Figure 2: 7^ for ensembles of surfaces obtained by blocking pure gravity triangula- 
tions k times. The values are plotted versus volume and the blocking increases in 
the left direction. 



by looking at their scaling behavior ^. The point-point distributions n{r, N) are 
defined as the number of nodes at a geodesic distance r from some random origin. 
Given that the system is close to criticality, it is plausible and indeed appears correct, 
that these distributions satisfy a scaling ansatz for large N. This implies that 

n{r,N) = Ny/ir/N^ = N'-'I''- f [-^) • (15) 

The scaling exponents x and y are fixed by assuming that there exist a power-law 
relation between a length scale I and volume N, oc l'^" , and by noting that the 
integral of n(r, A^) yields the total volume. 

In numerical simulations dn can be extracted by measuring n(r, A^) at dif- 
ferent volumes and find the optimal value that collapses all the data onto a single 
curve using Eq. (p!5|). The scaling assumption also implies that both the position tq 
and height n(ro) of the peak in n (r, N) should scale simply with volume N, i.e. 

r^^N'^H and n(ro) ~ "^h . (16) 

In Table || we show dn measured using the relations Eq. (|l6|) on volumes 
ranging from 250 to 1000 nodes. This is shown both for the initial (unblocked) 
lattices, and for ensembles obtained by blocking once or twice. It is clear that the 
extracted Hausdorff dimension for the blocked ensembles is completely in agreement 
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dninmax) 





3.15(2) 


3.13(1) 


1 


3.16(2) 


3.15(1) 


2 


3.15(2) 


3.20(1) 



Table 1: The Hausdorff dimension dn extracted both from the scahng of the location 
of the peak in n{r, N) and its maximal value. Results are shown for different number 
of blocking steps k with bi = 2. The volumes used are from 250 to 1000 nodes. 

with its value extracted from the bare lattice dn = 3.15(2). But note that these 
results are far from the exact result dn = 4. This effect has been observed before 
and is attributed to the presence of large corrections to scaling in the current lattice 
ensemble in which degenerate links are excluded. In Q it was shown that this 
estimate moves up to be statistically consistent with four if a wider ensemble is 
used in which tadpole and self-energy insertions are included. 

It is more instructive to compare the distributions directly. In Fig. ^ we plot 
distributions at volume N = 250, both for blocked and unblocked triangulations. 
This is compared to an estimate of the infinite volume distribution ?7-(r, oo). (the 
distribution n(r, oo) is obtained by scaling a distribution corresponding to a large 
volume {N = 2000) down to = 250 using Eq. ( p!5| ) with the exact value dn = 4). 
The finite size effects are very apparent in Fig. ^ as both distributions deviate from 
n{r, oo). But it is intriguing that the curve corresponding to the blocked = 250 
ensemble is closer to the estimated continuum curve - perhaps the first indication 
that the blocking is resulting in a flow towards a fixed point. 

Combined the measurements of 7^ and dn for the ensemble of blocked man- 
ifolds provide strong evidence that the fractal structure is preserved. This gives us 
confldence that the DB method correctly captures the long distance physics of the 
model. 

4.2 Fixed point behavior 

It is not enough that the RG transformation preserves the fractal structure, it should 
also induce a flow in the couplings towards a critical (unstable) fixed point. Opti- 
mally we would like to study the behavior of relevant operators around that fixed 
point and extract some critical exponents. Unfortunately for pure gravity the only 
simple relevant operator corresponds to the area operator, which is not accessible 
to us as we are working in the fixed area (volume) ensemble. 

Instead we have looked at the behavior of expectation values of various ir- 
relevant operators constructed out of the local coordination number qi or curvature 
(6 — Qi). Specifically we have looked at the following operators (where the super- 
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Figure 3: The point-point distributions n(r, A^) for = 250 in the case of pure 
gravity. Data for both blocked (87.5%) and unblocked triangulations is compared 
to an estimate of the continuum distribution n{r, oo). 



script k labels the blocking level): 



\ {ij) 



(18) 



(ij) implies that the sum is over neighboring nodes. 

If a fixed point exists then we would expect that in the absence of finite size 
effects (i.e. the initial lattice is lying on the critical surface) that these operators flow 
under repeated blocking into the fixed point. Thus differences of the expectation 
values between blocking levels should vanish; 



, /c ^ oo. (19) 



Since we are forced to work with finite volumes we have to be careful in how we 
compare these operator expectation values. It only makes sense to compare operator 
expectation values at the same volume, but differing in the number of blocking steps. 
This we do in Fig. ^. We have simulated at different initial volumes N^^^ (500 to 
8000 nodes) which then are blocked down to fixed volume (256 or 500 nodes) using 
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Figure 4: Differences in expectation values of the operator 0\ for surfaces obtained 
with different amount of blocking (but with the same volume) . 



14 

k 



a blocking factor hi = 5/4. Then the difference Eq. 19 is constructed. Clearly after 
just a few blocking steps this operator difference is close to zero. Similar results are 
observed in other channels and give strong evidence of the existence of a fixed point 
in the model. Combined with the previous demonstrations regarding 7^ and dfj it 
appears that the decimation method is a rather promising approach to a real-space 
RG for pure gravity. 



4.3 Higher derivative curvature term 

How robust is the ND method? To answer that we have studied the RG flow of other 
operators which are thought to be irrelevant for pure gravity. One such operator is 
O = log qi which corresponds to a set of higher derivative curvature interactions. 
The non-perturbative flows of this operator have been examined previously in the 
context of the LGB method [|l^]. The idea is to add such a term to the bare 
action with some coupling a. Under blocking the coupling to this new term should 
be driven to zero in the effective action as it corresponds to an irrelevant operator. 
Thus expectation values of operators should approach their values in the model 
without this extra interaction term. This is illustrated in Fig. ^ were we plot the 
RG trajectory in the (Oi,02) plane for a range of bare couplings a. 

We see that while the flows for different bare a start out from different 
points in the plane, for sufficiently positive a they are all attracted towards a fixed 
point. After passing close to this point the flows converge on a unique 'renormal- 
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Figure 5: The RG flow of the curvature operators Oi and O2 as an ensemble of 
triangulations is blocked. The flow is shown for different values of the bare coupling 
a. The dashed circle indicates an unstable fixed point, the solid a zero volume sink. 



ized' trajectory which ends at a trivial fixed point corresponding to zero volume 
triangulations. 

However, there are some indications that this no longer is true for sufficiently 
large negative values of a (the solid points). There the expectation values start to 
deviate considerable from the pure gravity values. At these lattice volumes they 
never pass close to the unstable fixed point but run smoothly to the trivial end 
point. Such an effect has been observed before and it is believed that for large 



negative a the model enters a 'super-crumpled' state |11|. 

Such an effect is also observed in the fractal structure. We have measured the 



a = 10 
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a 



-1 



a 



a 



-10 



1000 

1 500 

2 250 

3 125 



-0.47(1) -0.47(2) -0.48(1) -0.49(1) -0.03(3) 2.1(1) 

-0.49(1) -0.47(2) -0.50(1) -0.48(1) -0.30(2) 1.8(2) 

-0.51(2) -0.46(2) -0.47(2) -0.48(2) -0.53(2) 2.0(2) 

-0.47(1) -0.42(2) -0.46(2) -0.42(1) -0.65(3) 0.5(2) 



(k) 

Table 2: 7s for ensembles of triangulations obtained by blocking k times. Results 
are shown for various values of the bare coupling a. 
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string susceptibility exponent as before for several values of a and various number of 
blocking steps. The result is shown in Table For 10 > a > —1 we get consistently 
the pure gravity value. But for a = —5 the initial ensemble yields a different value 
of 7s; but as it is blocked the value appears to flow towards the pure gravity fixed 
point value. For a = — 10 the fractal structure is dramatically changed, we get a 
consistent value of 7^ ~ 2, even after considerable blocking. 

It is not clear at present whether these results can be explained on the basis of 
finite size effects which are much more dramatic for large negative a or whether they 
indicate that no continuum limit exists for such models - the trajectories are never 
attracted towards the usual unstable fixed point. Certainly we have not observed 
any non-trivial new fixed points associated with these negative coupling flows. 

5 Ising model coupled to 2d gravity 

The next step is to apply the method to matter fields living on the surfaces. For that 
purpose we chose Ising spins coupled to gravity in which case the matter contribution 
to the partition function Eq. (|lO|) is 

Zm{T) = Y^eP^{ij)'^i'^3 . (20) 

Here the sum runs over each edge length in a triangulation and the coupling (3 is 
tuned to its critical value. 

This model is exactly solvable [^] so we know the critical coupling, /3c = 
^ log together with the exponents governing the critical behavior. In addition we 
know the effect the matter has on the fractal structure. The Ising spins only change 
the string susceptibility exponent at the critical point to 7^ = —1/3, elsewhere it 
retains the pure gravity value. 

We have simulated this model at the critical point of the Ising model for 
lattice volumes = 500, 1000 and 2000 nodes. A Swendsen-Wang cluster algorithm 
was used to update the spins and again about 10*^ sweeps performed per lattice 
volume. 

5.1 Blocking the spin configurations 

When it comes to extending the ND method to block a spin configuration, in ad- 
dition to a triangulation, we have to decide to what extent the blocking of the 
gravitational sector and spin sector should influence each other. The simplest so- 
lution is that the blocking of the triangulations should be performed precisely as 

^The actual value of 7^ should not be taken too seriously as the baby universe distribution for 
a = — 10 does not fit very well the asymptotic form Eq. (p^. But for a crumpled phase there is 
no reason to expect such a simple asymptotic behavior, as has been observed in simulations of 4d 
simplicial gravity 
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for pure gravity - the only effect of the critical matter will be to alter the initial 
distribution of triangulations. When a given blocking ratio has been obtained a new 
spin configuration is created from the spin configuration of the initial triangulationQ. 

To block the spin configuration we choose the "majority-rule" transformation 
that has been applied successfully to the Ising model on a regular lattice. The spins 
on neighboring nodes are grouped into "blocks" and a block spin value assigned to 
each block depending on the sum of its spins. If the sum is positive the value is +1 
but -1 if it is negative. If the sum is zero we assign +1 and -1 with equal probability. 
However, there is a difference in applying this method to random surfaces as opposed 
to a regular lattice. A regular lattice can be divided into equal non-intersecting 
blocks, an arbitrary triangulation not. Instead we define a "block" as a spin that 
survives the decimation together with the neighbors it had before blocking. This 
definition is appealing as it incorporates the influences of the geometry. The size of 
the "spin-blocks" will be depend on the local curvature. 

However, there is still an ambiguity in this definition of "spin-blocks". Some 
spins will contribute to more than one such block, in particular if two neighboring 
nodes both survive the decimation. This is especially true if a small blocking ratio 
is used. Thus it might be necessary to limit the influence of such nodes by excluding 
them from other blocks than their own. We tried this modified definition of spin- 
blocks, but it did not change the results appreciably. 

5.2 The fractal structure 

How well does the node decimation preserve the fractal structure for this model? 
Again we have computed the string susceptibility exponent; the result is shown in 
Fig. ^ for three different initial volumes. Notice that as the blocking of the geometry 
is independent of the spins the results do not depend on the definition of spin-blocks. 

The first thing to notice is that without any blocking we measure the value 
-fg ^ —0.41, which is somewhat lower than the theoretical value 7s = —1/3. This 
is due to the fact that we measure at the infinite volume critical coupling not 
the pseudo-critical coupling appropriate for the finite volume used. This effect 
has been observed in earlier numerical estimates of 7^ for this model in ||^]. There a 
sharp peak was observed around j3Y with peak value 7^ —1/3. But away from the 
peak the values rapidly return to the pure gravity value 7s = —1/2. As the volume 
is increased the peak moves closer to the true critical point. 

However under blocking the measured 7^ appears to flow towards the exact 
value. This, of course, is very encouraging from the point of view of capturing fixed 
point physics. Notice that direct simulations of the smallest volume lattices would 
have produced estimates for 7s further from the exact result. 

We tried methods in which the surface and the spins are blocked simultaneously. After re- 
moving each node a new spin configuration is found and in addition the possibility for choosing a 
given re-triangulation is allowed to depend on the spins. These methods did not seem to capture 
the physics of the system correctly. 
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Figure 6: Changes in 7^ as triangulations with Ising spins are blocked. Data is 
shown for three initial volumes. Again the number of blocking steps k increases to 
the left. The horizontal line is the exact value 7^ = —1/3. 



5.3 Scaling exponents 

To measure the critical exponents of the Ising model wc consider several operators 
constructed out of products of neighboring spins and apply the methods described 
in section 2. In practice the possible operators can be classified according to their 
symmetry properties under the Z{2) operation cjj —cTi- In the linearized matrix 
Tq,^ the cross-correlations between these even and odd operators vanish. 
In the odd sector wc have looked at the operators 

SI = and = ^(7^(7^(7^, (21) 

i (ijk) 

where the first sum is over nodes and the second over triangles {ijk). For the even 
sector we have chosen the operators 

Sf = ^aiaj , SI = ^GkCTi and SI = aiGjOkai. (22) 

{ij) {kiy {ijkl) 

The indices i,j, k, I refer to an elementary link ij together with the complementary 
vertices k, I associated with the two neighboring triangles. 

These operators allow us to determine numerically two scaling exponents; 
the magnetic exponent yi and the thermal exponent y2- The former corresponds to 
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the odd operators, the latter to the even. These scaUng exponents are related to 
the critical exponents of the Ising model v and b through the relations 



y\ = T— r and 2/2 = (23) 

+ 1 UClH 

For the Ising model coupled to gravity these critical exponents are known, i/dn = 3 
and 6 = 5, yielding the values yi = 5/6 and 7/2 = 1/3 for the scaling exponents. 
In comparison, for the Ising model on a flat regular lattice these exponents are 
yi = 15/16 and y2 = 1/2. 

In Figs. 0a and b we show the extracted values of the scaling exponents as a 
function of the RG iterations k using a blocking ratio bi = 2 (and for different initial 
volumes). The values shown are obtained including only one spin operator for each 
sector in the calculation of the matrix Tq,^, i.e. 3° and Sf. With the accuracy we 
have at present no obvious improvement was observed by enlarging the operator 
basis (although the results were perfectly compatible). 

The results show that after one or two blocking steps we have approached a 
critical fixed point. It clearly is the correct one as we get results which agree with 
the exact values within 2 or 3% accuracy. In comparison the exponents for the Ising 
model on a flat lattice are clearly ruled out. It is impressive that even after blocking 
down to triangulations consisting only of 16 nodes the critical exponents of the Ising 
model still come out correctly. 

Do these results depend on the blocking ratio? Unlike for the triangulations 
the blocking of the spin configurations is not independent of the amount of node 
decimation. We tried three different blocking ratios, bi = y/2, 2 and 4, but the results 
did not change very much. The error-bars increased somewhat as the blocking ratio 
was decreased, which is not unexpected. 



6 Alternative blocking methods 

As mentioned in the introduction there are other real-space RG transformations that 
have been proposed for random surfaces. But how well do those methods capture 
the long distance physics of the surfaces? We have investigated this for the local 
geodesic blocking (LGB) method. 

A brief description of the method given in |^]: Given a triangulation T{N) we 
chose a subset of N' nodes and connect them together to form a block triangulation 
T'. This is done in such way as to preserve locally the relative geodesic distances 
between block nodes. In practice the subset of nodes is fixed at the beginning of the 
simulations. The initial triangulation is then evolved using a standard link flip algo- 
rithm and simultaneously the block triangulation is updated with the prescription 
that a link lij is flipped to Iki if the geodesic distance dki is less than dij (measured 
on the unblocked lattice). That this constraint is local is essential for the algorithm 
to be practical in numerical simulations. 
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Figure 7: The scaling exponents yi and y2, corresponding to odd and even spin 
operators respectively, versus the number of RG iterations k. In both plots the 
lower horizontal lines are the exact values for an Ising model coupled to gravity 
yi = 5/6 and y2 = 1/3, while the upper lines are the flat space values yi = 15/16 
and 2/2 = 1/2. 
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bi 


k = 

N 7s 


k = 1 

N 7s 


k = 2 

N 7s 


Sphere 

1.08 
1.2 
1.6 
2.0 
4.0 


270 -0.504(14) 
300 -0.506(13) 
400 -0.486(16) 
500 -0.484(17) 
1000 -0.499(11) 


250 0.271(18) 
250 0.280(11) 
250 0.296(27) 
250 0.273(30) 
250 0.289(40) 




Torus 

4 
4 


2304 -0.486(16) 
1024 -0.504(9) 


576 0.306(27) 
256 0.268(35) 


144 0.281(30) 



Table 3: Measured values of 7^ for triangulations blocked with the LGB method. 
Results are shown for both spherical and toroidal topologies and for various blocking 
ratios 6;. 



We have measured the distribution of baby universes on the ensemble of 
triangulations created by this blocking method, both for surfaces with spherical and 
toroidal topology^ (the latter was used in Q). The result for 7^ are shown in Table 
3. They indicate that the fractal structure is not correctly preserved. The ensemble 
of blocked triangulations has, consistently, 7^ ~ 0.28 independent of blocking ratio 
or topology. 

This is in our opinion a serious flaw of the method. If 7^ is thought of as a 
regular critical exponent this result would imply that this RG scheme is driving the 
system towards a new fixed point - not pure two-dimensional gravity. A possible 
explanation is that as the method uses local constraints to dictate the link flips on 
the block triangulations, it may actually simulate the effect of having some kind of 
matter field living on the manifolds. This could effectively change the physics of the 
model. 

The situation might be improved by using a generalization of the LGB 
method; a global geodesic blocking (GGB). It is defined by requiring that the 
block triangulation preserves (up to an overall global length rescaling) as closely 
as possible all possible geodesic paths between block nodes. In practice this may be 
implemented by minimizing a cost function E when performing block lattice link 
fiips 

E = J2{d,,-Xd!,^y . (24) 

^Measuring 7^ for manifolds of toroidal topology is slightly more complicated than for a sphere. 
Firstly, the minimal loop used to identify a baby universe must not be topologically non-trivial. 
Secondly, the baby universe distributions now describe spheres growing on a torus. Hence there are 
two different -fs exponents in Eq. (jl^), one for each topology. For pure gravity both these exponents 
come out correctly in simulations. 
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The quantities dij and d'^j are the distances in units of the cut-ofF between block 
nodes i and j defined on both the initial lattice and block lattice respectively. The 
parameter A relates the two scales and may be chosen so as to minimize E for 
fixed lattices. Clearly other choices of E are possible corresponding, for example, to 
weighting small and long geodesies differently. 

We have investigated this method numerically, using the cost function Eq. 
(24) and found that it indeed seems to cure the problem of an incorrect 7^ exponent. 
The values of 7^, obtained with limited statistics, are consistent with the pure 
gravity value. Unfortunately as the method is highly nonlocal it is computationally 
burdensome and not feasible in practice. 



7 Outlook 

We have presented results on a new real-space renormalization group method for 
dynamically triangulated 2d gravity. Our results are encouraging in so far as they 
indicate that under a simple node decimation (ND) flows toward a fixed point in 
the critical surface are produced. Measurements of the string susceptibility and 
Hausdorff dimension strongly support the notion that this fixed point corresponds 
to 2d quantum gravity. The effect of including an irrelevant operator corresponding 
to a higher order curvature term also support this picture. Indeed the demonstrated 
irrelevance of such an operator outside of perturbation theory might be taken as the 
first successful non-trivial application of the method. 

We have extended the method to a simple critical Ising system coupled to 
2d gravity and have been able to extract exponents in excellent agreement with 
the predictions for gravitational dressing given by the DDK/KPZ formula |14|. Fur- 



thermore measurements of the string susceptibility exponents indicate that the finite 
size effects are significantly reduced. This provides further evidence for a flow to 
the correct fixed point. 

There is still some scope for improving the method. A more systematic 
investigation of how to inter-disperse the node and spin blocking would be desirable. 
Also the effects of employing a more extensive operator basis should be studied. 

Provided the method or refinements thereof survive these tests the way is 
open to apply these RG methods to a variety of unsolved problems in 2d gravity. 
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